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Abstract. In this paper we study phase transitions for the Potts model with 
four competing interactions (external field, the nearest neighbor, the second 
neighbors and triples of neighbors) on Cayley tree of order two. We prove 
that for some parameter values of the model there is phase transition. We 
reduce the problem of describing limiting Gibbs measures to the problem of 
solving a system of nonlinear functional equations. We extend to the Potts 
model setting the results obtained by Ganikhodjaev and Rozikov %_ on phase 
transition for the Ising model. 
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1. Introduction 

It is well known that lattice spin system is a large class of systems considered 
in statistical mechanic. Also it is well known that the structure of the lattice plays 
an important role in examinations of spin systems. The Potts model, which was 
introduced by Potts [in 1952] as a generalization of the Ising model, has received an 
increasing theoretical and experimental aspects in recent years (see [T3] for details) . 
The idea came from the representation of the Ising model as interactions spin which 
can be either parallel or anti parallel that is up or down. 

We consider a semi-infinite Cayley tree J*^ of order k >2, i.e., a graph without 
cycles with (k-l-1) edges issuing from each vertex except a;° and with k edges issuing 
from the vertex x^, which is called the tree root. We generalize the problem in 
[3] for the Potts model with four competing interactions. This problem has been 
investigated for the Potts model with two competing interactions (see ^ for details). 

There are many approaches to derive equation or system equations describing 
limiting Gibbs measure for the lattice models on Cayley tree (see [I], [10], [TT] . 
[12]) [13] for details). One of these approaches is based on recursive equations for 
partition functions ([13]). Another approach is Markov random fields (see [4], [8]). 
Naturally both approaches lead to the same equation [J. We concentrate on the 
first approach, that is, we study recursive equations for partition functions. 

In this paper we study the problem of phase transition for the Potts model 
with four competing interactions (external field, the nearest neighbors, the second 
neighbors and triples of neighbors) on Cayley tree of order two. In [S], it has been 
proved that there is phase transition for some parameter values of the Ising model. 
We generalize this problem for the Potts model. To construct the Hamiltonian 
equation we consider the generalized Kroneker's symbol S well coordinated with 
theory of quadratic stochastic operators. 
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The rest of the paper is organized as following. In section 2 we give basic defi- 
nitions and notations. In section 3 we obtain the recursive equations for partition 
functions. Section 4 contains the solutions of the system of nonlinear equations. 
The solutions describe whether there is phase transition. 

2. Preliminaries 

Cayley tree. The Cayley tree F'' (see [T]) of order fc > 1 is an infinite tree, i.e., 
a graph without cycles, from each vertex of which exactly fc + 1 edges issue. Let 
r*^ ~ (V, L, i), where V is the set of vertices of F*^, L is the set of edges of F*^ and i 
is the incidence function associating each edge £ € L with its ends points x,y €V. 
If = {x,y}, then x and y are called the nearest neighboring vertices and we 
write i =< x,y >. For x, y €z V, the distance d{x, y) on Cayley tree is defined by 
formula 

d{x, y) = min{d\x = xq, xi,X2, Xd-i, Xd — y (z V such that the pairs 
< Xo,Xi >, < Xd-i,Xd > are neighboring vertices}. 
For the fixed x'^ £ V we have 

Wn = {x£V : d{x°,x) = n} 
Vn = {xeV : d{x°,x) < n}, 
Ln ^ {£ ^< x,y >e L : x,y e 14}, 
A collection of the pairs < xo,xi >, < Xd-i,y > is called a path from x to y. 
We write x < y ii the path from x° to y goes through x. We call the vertex y a 
direct successor of a;, if y > a; and x, y are the nearest neighbors. The set of the 
direct successors of x is denoted by S{x), i.e.. 

Six) ={y€ Wn+i : d{x, y) = l,xe W„}. 

It is clear that for any vertex x 7^ x", x has k direct successors and has k + 1. 
The vertices x and y are called the second neighbor which is denoted by > x, y <, 
if there exists a vertex z £ V such that x, z and y, z are the nearest neighbors. 

The second neighbors > x, y < is called one-level neighbors, if vertices x and y 
belong to Wn for some n, that is, if they are situated on the same level. Three 
vertices x,y and z are called a triple of neighbors, which is denoted by {x,y,z), 
(x,y), {y,z) are nearest neighbors and x ^ z. We will consider only one-level the 
second neighbors. 

We assume for neighbors triple (x, y, z). For neighbors triple (x, y, z) we consider 
the following Kroneker's symbol 

{1, if a{x) = a{y) = a{z) 
i, if a{x) ^ a{y) ^ a{z) or a{x) ^ a{y) = cr{z) 
0, else 

where, x,z £ Wn for some n and y G Wn-i- We can write (2.1) as follows: 

This definition is well coordinated with theory of quadratic stochastic opera- 
tors, where quadratic stochastic operator corresponding to generalized Kroneker's 
symbol in (2.1) is identity transformation [9|. 
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In the Potts model, the spin variables (j{x) take their values on a discrete set 
$ = {1,2, ...,(7}, q > 2 which are associated with each vertex of the tree (see [6] 
for details). 

In this paper, we consider the Potts model with four competing interactions is 
defined by the following Hamiltonian 

{x,y) >x,y< 

(2.2) — J2 ^ Sc(x)cr{y)cr(z) ~ h ^ Soa(x) 

{x,y,z) x£V 

where the sum in the first term ranges all nearest neighbors, second sum ranges 
all the second neighbors, the third sum ranges all triples of neighbors and the spin 
variables (t{x) assume the values 1, 2, 3. S is the Kroneker's symbol and J, Ji, 
J2 and h E M. are constant. We assume that J.Ji.J2-h ^ 0. In |6J, we solved the 
problem of the phase transition described by means of the equation (2.2) for h = 
and Jl = 0, it has been found exact solution for the problem of phase transition 
in the Potts model with competing ternary and binary interactions with external 
field. 



3. Recursive equations for Partition Functions 

Let A be a finite subset of V. We will denote by 17(A) the restriction of a to 
A. Let ct(V\A) be a fixed boundary configuration. The total energy of cr(A) under 
condition ct(V\A) is defined as 

^a(ct(A)|ct(1/\A)) = - J 5a(x)a{y) - J ^ ^cr{x)a{y) 

<x,y>:x,y£A <x,y>:xeA,y(A 

>x,y<:x,y£A >x,y<:x£A,yfA 

(3.1) - J2 Y ^aix)aiy)a{z) - J2 ^ <y(x)a(y)a(z) 

<x.y,z>:x,y,z£A <x,y,z>:x,z^A,y£A 

— h Y^ SQcr(x) 

xeA 

Then partition function Zf^{a{V \ A) in volume A under boundary condition 
a{y\h) is defined as 



(3.2) Z^{a{V\K))^ Y exp(-/3i/A(a(A)|a(nA))), 

cr(A)en(A) 

where r2(A) is the set of all configuration in volume A and /3 = ^ is the inverse 
temperature. 

We consider the configurations cr(V„) and the partition functions Zy^va. volume 
V„ and for the brevity it is denoted them as (t„ and Z^"^ respectively ([4], [6]). Let 
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US decompose the partition function Z^"^ into following summands 
where 

(3.3) Z\^^ = eM-PHvMnV\V.n))). 

Let 9 = exp{pj), 9i — exp{pji), 02 = exp{pj2) and 6*3 = exp{j3h). 

In this paper, we consider q — i and fc = 2. It should be noted that in [3] we 
have ignored the all triples of neighbors. 

Let S{x°) = {a;\ x^}. If a{x^) = i, cr(x^) = j and cr(x2) = m, then from (|3?T|) 
and (|3.2p . Thus, we have following equations; 
3 

j,m=l 

so that 



^3 



((Zf + (^3 ) ) + '^zl^'^'zl 
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and 

After replacement u„ = ^7;^ and Vn ~ —j:^ we have the following system of recur- 
rent equations 

0, + 20^Un-l{l + Vn-l) + 2W„_1 + 0^9i92ul_^ + 0ivl_^ 



9i02 + 20^2{Un-\ + Vn-l) + 6*1 ("Ll + ^l-x) + 2m„_iW„-1 



0-^0^02 + 20V02{Un-l + Vn-l) + + f^l) + 2u 

If u = limUn and u = limvn then we have the following; 

6li + 20yW2u{l +V) + 2V + 0^0i02U^ + 0iv^ 

03U = 

(3.4) 03V = 



0^0102 + 20V02{U + v)+ 0i{u^ + w2) + 2UV 
0l + 2u + 20y%vil +U) + 0iU^ + 0^0i02U^ 
0^0102 + 20V02iu + v)+ 6*1 (m2 + w2) + 2uv 



The number of the solutions of equation (|3.4p depends on the parameter (3 — -j^- 
The phase transition usually occurs for low temperature. It is possible to find an 



PHASE TRANSITIONS FOR THE POTTS MODEL WITH FOUR COMPETING INTERACTIONS5 



exact value of T* such that a phase transition occurs for all T < T* where T* is 
called a critical value of temperature. 

Finding the exact value of the critical temperature for some models means to 
exactly solve the models. 

4. Translation-invariant Gibbs measures 
In this section we describe the solutions of the system of nonlinear equations 
9i + 20y/9^u{l + v) + 2v + 6^eie2U^ + Oiv^ 



nu = 



(4.1) e^v = 



e^ei02 + 2eV02{u + v) + 0i{u^ + v^) + 2uv 

6i + 2u + 2B^f¥2v(\ + m) + Qxu^ + d'^dxQ2V? 



9^6162 + 2eV02{u + v) + 9i{u^ + u2) + 2uv 

Note that if there is more than one positive solution for equations ()4.1|) . then 
there is more than one translation-invariant Gibbs measure corresponding to these 
solutions. If the equations ()4.ip have more than one positive solution, then it is 
said that a phase transition occurs for the model (|3.2p ([5], [5]). 

Let us solve the equations (|4.ip . Subtracting these equations we have 

(4.2) ^^^^ _ _ 2(0 - l)iu-.)+ 0,(0^ -Diu^-.^) 



+ 29{u + v) + 2uv + 6*1(712 + t;2) 



so that equality u = v give us some solutions of (|4.ip . where 9 = 9y/92. For 
u — V we have 

,^ , „ {9^9i+29 + 9i)u' + 2(9 + \)u + ei 

(4.d) 9^u — = = 

2{9i + l)u^ +Wu + 9^9i 

Let us consider function 

, {9''9i+29 + 9i)u^ + 2{9 + l)u + 9i 

(4.4) Jiu) = = = . 

2(6ii -f l)u2 + 46)u + 6»26li 

In order to describe the solution of the system of nonlinear equations in the (|4.ip . 
we should analyze the function in the (|4.4p . 

Now we can study elementary analysis of this function as following: 

1) Let us take the first derivative of /(u), in this case, if 6* > 1, that is 92 > 
then f{u) is increasing and if < 1, f{u) is decreasing for u > 0; 

2) Let us take the second derivative of f{u). If 9i > 91, where 



_9^ + 9 + l + V99* + 2603 + 35^2 + 509 + 33 
' " (^2 + 2)(^+l) ' 

then exist inflection point u* > such that /"(m) > for < u < u* and f"{u) < 
for u > u* . 

Thus we can generalize the Lemma proved in ^31 as follows; 

Lemma. Let 9 > 1, that is 6*2 > ^ and there is inflection point. There exist 
771(6', 6*1, 6*2), r]2i9, 6*1, 6*2) with < 771(6*, 6*1, 6*2) < 772(6*, 6*1, 6*2) such that the equation 
(|4.ip has three solutions if 771(0, 6*1, 6*2) < 9^ < ^72(0, ^i, 6*2); has two solutions if 
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either 63 — i]i{9, 61,92) or 63 = ??2(6', ^1, 6*2) and has single solution in other cases. 
In fact 

where ul,U2 are the solutions of equation 

uf (u) = f{u). 

Proof follows from properties 1) and 2) of function /. Equation uf'{u) — f{u) is 
equivalent to following equation 

(4.5) Au^ + Bv!^ + Cu^ + Du + E = 0, 

where A = 2(0i + l){6^ei + 20 + 0i) > 0, B = 8(0 + l)(0i + 1) > 0, Z? = 800i > 0, 
E = 0201 > 0. 

As C = -e^el - 20103 - 0202 + 802 + 80 + 60? + 601, 

then equation ()4.5p has no positive roots if C > 0, and have two positive roots 
1*2 if C < 0. It is easy to show that exist 0" such that for all 0i > 0** coefficient 

C is negative, where > \ -g-- Then we have the following Corollary. 



Corollary. For 0(02) > , 0i(0,02) > max{0l, 0^) and r;i(0, 0i,02) < 

^'3 < &i,d2) the equation (|4.5p has three positive roots. Thus there is phase 

transition. 

Now let us come back to equation (|4.2[) . After cancelling to (u-v) we have 

^ 0102^2 + (2020 ^ 01 (02 - 1)).S + 020102 - 2(0 - 1) 

^'■'^ ' = ' 

where u + v — s and uv = t. After dividing the first equation of (|4.ip to the second 
and simplification wc have 

(4.7) ^^0i£l+2^ 

0201 + 01 - 20 

From (|4.6|) and (|4.7|) we obtain the following equality; 



03 0i(02 + l)-20 

01 [01 (03 + 02 + 20 + 2) - 2(02 + + 1)] 
0i(02 + l)-20 

It is easy to show that if 81(0,62) > then all coefficients in (|4.8p are 

positive. 

Note that the system of equations 



u + v = s 

(4.9) uv = t 

2 

have solutions when i < so that from (|4.7p follows 

(4.10) (0i(02 - 3)- 20)s2-8s-40i >0 



0i(0 + l)g + 2 _ 0i[0i(0 + l)-2]s2 + 2[0i(02 + 0-2)-2(0 + l)]s 

03 

(4.8) 
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Then there is 



4 + ^/16 + 46ii(6ii(6i2 - 3) - 261) 

> 0, 



6i6^ - Wi - 20 
such that for all s > s* is valid. 

Now let us consider equation (|4.8p . The elementary analysis get following: 

1) If 

03 < 

then equation (|4.8p has 

2) If 

O3 > 



2 [01 (^2 4 


-1) 


- 29] 




0i[ei{9^ + 9"^ + 29 ~ 


^2) 


-2{9^ + ( 


^+1)]' 


single root. 








2[9i{9'^- 




- 29] 





9i[0i{9^ + 9^ + 29 + 2)- 2(6*2 + e + l)] 
there is /cg such that the line V~ — feps is a tangent for parabola from right-hand 

side of (|4.8p and then for 63 < the equation (|4.8p has two positive roots 

and only one of them greater than s* . 

Thus the equation ()4.8|) has one positive root greater than s*, when 



> 

and 

2[9,{9^ + l)-29] ^ gi(g'+l) 

0i[0i(^3+^2 + 2^+2)-2(^2 + ^+i)] ^ fco ' 

where /cq is defined above. 

Recalling results for equation (|4.3p . we have the following: 

the system of the equation (|4.ip has five solutions if 6* > ^/j-, 6*1 > A and £? < 



6*3 < C, where 

A = max(6li,6li*) 

2[6'i(^2 + 1) - 2^1 

(4.11) S = max(7yi,- L iv "r ; j 



?i [9i (613 + 612 + 26» + 2) - 2(02 _^ _^ 1)] ' 



C = mm(ry2, ; ). 

fco 



fco - 

More detail analysis similar [2] and [3] get that only three solutions are stable. 
Thus we have proved the following theorem. 

Theorem. Assume the conditions (I4.1ip are satisfied then for the model (13. ip 
there are three translation- invariant Gibbs measure, i.e., there is phase transition. 



It is easy to show that these three measures correspond to boundary conditions 
a{V\V,,) = i, i = l,2,3. 
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Remark. We can similarly study the solutions of new system of nonlinear 
equations by using Kroneker's symbol 6 that has the form 
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